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Abstract: We construct an algorithm to determine all stationary axi-symmetric 
solutions of 3-dimensional Einstein gravity with a minimally coupled self-interacting 
scalar field. We holographically renormalize the theory and evaluate then the on- 
shell action as well as the stress tensor and scalar one-point functions. We study 
thermodynamics, derive two universal formulas for the entropy and prove that global 
AdS provides a lower bound for the mass of certain solitons. Several examples are 
given in detail, including the first instance of locally asymptotically flat hairy black 
holes and novel asymptotically AdS solutions with non-Brown-Henneaux behavior. 
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1. Introduction 

Exact solutions have a long history in General Relativity, going back all the way to 
Schwarzschild in 1916, see [1] and references therein. Of particular interest are axi- 
symmetric stationary solutions: they have two commuting Killing vectors, both of 
which are phenomenologically viable symmetries in many physical situations, ranging 
from applications in astrophysics (black holes) to applications in the gauge/gravity 
correspondence (reasonable gravity duals for ground states or meta-stable states of 
the dual field theory). 

Given some gravitational theory with matter — for the sake of concreteness let 
us restrict to Einstein gravity — it is notoriously difficult to derive all solutions to 
the equations of motion (EOM), which are non-linear coupled partial differential 
equations (PDEs). In 3+1 dimensions the imposition of axi-symmetry and station- 
arity still leads to quite complicated equations of that type. In 2+1 dimensions the 
PDEs are converted into ordinary ones, which is an essential simplification; however, 
even this simplification, which we implement in the present work, is in general not 
sufficient to permit an efficient algorithm to find all solutions. 

It is much simpler, if not trivial, to find solutions in an inverse way: take some 
metric, declare it to be a solution of the Einstein equations with matter and deduce 
what kind of matter is needed for an exact solution. However, this procedure typically 
leads to solutions devoid of physical interest. 

In this paper we propose a way to find all stationary axi-symmetric solutions to 
2+1 dimensional Einstein gravity with a scalar field that simultaneously addresses 
the issue of physical relevance, similar in spirit to "designer gravity" [2]. 

Our main ingredient is to avoid fixing the scalar self-interaction potential, but to 
extract it as an output. The crucial input, instead, is a certain function that deter- 
mines essential properties of the geometry: its asymptotic behavior, the number and 
type of Killing horizons, the singularities etc. Thus, one can design the desired kind 
of geometry, while simultaneously guaranteeing the presence of reasonable matter — 
a minimally coupled self-interacting scalar field that is regular, at least outside a pos- 
sible black hole horizon. We call this procedure "bootstrapping gravity solutions". 
As we shall demonstrate, our procedure is particularly useful in the context of the 
Anti-deSitter/Conformal Field Theory (AdS/CFT) correspondence [3]. Let us state 
now some of our main results. 

• One key result is our solution generating algorithm presented in section 2.1, 
which requires one input function /i subject to an inequality, and several inte- 
gration constants. Appropriate choices for the input function and integration 
constants are severely constrained by physics requirements (regularity, smooth- 
ness, suitable asymptotic behavior, absence of closed time-like curves, ...) dis- 
cussed in section 2.2. 
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• We find specific 2-parameter families of stationary black hole solutions in sec- 
tion 2.3 and show that each such family has exactly one associated regular static 
solitonic solution. Moreover, we identify a constant of motion that generically 
leads to singular solutions when it is varied. 

• In section 3 we construct explicitly a family of locally asymptotically flat black 
hole solutions with scalar hair, with a symmetry-breaking scalar potential de- 
picted in figure 1. We also prove that there are no globally asymptotically flat 
black hole solutions with scalar hair. 

• For asymptotically AdS boundary conditions (relaxing Brown-Henneaux) we 
holographically renormalize the theory in section 4. We calculate the on-shell 
action and 1-point functions for Dirichlet, Neumann and Mixed boundary con- 
ditions for the scalar field. We focus particularly on marginal multitrace defor- 
mations, since in this case the 2-parameter black hole families that we construct 
can be interpreted as states of the same boundary theory. 

• We generalize the results of section 4 and study the most general asymptotic 
AdS conditions compatible with locality of the holographically dual CFT in 
section 5. We specifically address subtleties that arise for resonant scalars and 
Breitenlohner-Freedman (BF) saturation. 

• We analyze the thermodynamics of black hole solutions in section 6 and prove 
two simple formulas for the entropy, one of which has a natural interpretation 
in terms of the Cardy formula, while the other one proves a 3-dimensional 
version of Penrose's Weyl curvature conjecture. Moreover, we discuss the ther- 
modynamic (in) stability of black hole solutions and their associated solitons 
depending on the boundary conditions. 

• In section 7 we show how to recover known solutions within our algorithm and 
provide several new ones. We present in closed form two families of Brown- 
Henneaux solitons, with a potential displayed in figure 2, as well as non-Brown- 
Henneaux solutions with oscillating potentials exhibited in figure 3. Since the 
most difficult step in our algorithm is the integration of a linear first order 
ordinary differential equation actually arbitrary solutions can be constructed, 
albeit in many cases one has to resort to (simple) numerics. In our set of 
examples we have cherry-picked those that can be presented in terms of fairly 
elementary functions. 

In section 8 we conclude with possible generalizations of our algorithm. 

Before starting we mention some of our conventions. We use signature (— , -|-, +) 
and fix Newton's constant by IGttGn = 1- 
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2. Einstein gravity with scalar field 



Consider Einstein gravity with a minimally coupled self-interacting scalar field, to 
which we shall refer as Einstein-dilaton gravity (EDG). The EDG bulk action is 



S 



J d'xy^ + 2 - ^ i9(t>)' - Vicj))] (2.1) 



where V{(f)) is some (unspecified) self- interaction potential for the scalar field 0. The 
potential is defined in such a way that V{(f)) = leads to a negative cosmological 
constant with unit AdS radius. 

We are interested in stationary axi-symmetric solutions of the field equations 
descending from the action (2.1). To this end we employ a method devised by 
Clement [4]. Namely, the line-element is parametrized as 



ds^ = X., dx" dx" + " (2.2) 



^..= ^1/^"^.) (2.3) 



with the matrix A 

Y T-X, 

and all coefficient functions T, X, Y as well as the Einbein e only depend on the 
coordinate p. It is convenient to introduce the Minkowskian target space vector 
X = (T, X, y), whose norm is determined as X^ = — det A. The parametrization 
(2.2) is accessible for any stationary axi-symmetric solution. Inserting it into the 
action (2.1) and dropping a boundary term yields the reduced action 



S = Vo\ / dpe 



X +2 X'(l)' -V{(j)) (2.4) 



PO 



where 'Vol' refers to the volume of the 2-dimensional spacetime that was integrated 
out, dots denote d/dp derivatives, and the integration domain is bounded by some 
p = Po and p = oo. 

Solutions of the reduced action (2.4) are in one-to-one correspondence to station- 
ary axi-symmetric solutions of the original action (2.1). After varying the reduced 
action with respect to the fields e,X,(f) we gauge-fix the Einbein to unity, e = 1. In 
this gauge the field equations read 

Einstein: X = -^^X (2.5a) 

Klein-Gordon: X^^ + 2XX<j) = V'{(f)) (2.5b) 

Hamilton: - ^ X% 2 + ^ X^^^ - 1/(0) = (2.5c) 
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A first integral of the EOM is provided by the conserved particle angular-momentum 
(not to be confused with the black hole angular momentum of specific solutions), 



eijkX^X'^ = J k = const. (2.6) 

The Klein-Gordon equation (2.5b) can be reformulated as an equation for V , which 
can then be integrated. The integration constant is fixed by the Hamilton constraint 
(2.5c), which then leads to the following representation of the potential as a function 
of p. 

1 d2 

4 dp2 



np) = 2-7:r;2^' (2-7) 



2.1 Solution generating algorithm 

The Einstein equation (2.5a) implies planarity of the vector X. Thus, with no loss 
of generality we parametrize generic solutions to the field equations (2.5) as 

X(p) = Xi/i(p) + X2/2(p) (2.8) 

with some constant vectors Xi,2- The action (2.4) is Lorentz invariant, a property 
inherited from the 5'L(2,M) invariance of the parameterization (2.2), which we can 
use to bring the vector X into some convenient form. In fact, it is always possible 
to find Lorentz transformations that simultaneously set to zero the F-component in 
both vectors Xi 2- This is seen most easily if at least one of the vectors X12 is time- 
like: Go to the corresponding vector's rest-frame by appropriate Lorentz boosts, so 
that the X- and F-components vanish; then rotate in the XF-plane to eliminate the 
F-component of the other vector. Similar constructions work for all other cases. ^ 
Thus, with no loss of generality we rewrite the parametrization (2.8) as 



X{p)= \x\=-\-h-fA (2.9) 





which implies that the 2-dimensional metric A^,^ in (2.3) simplifies to 

/i 



^Onc can show this for the cases where the vectors are space-hke or hght-hke using successively 
Lorentz transformations. For instance, we can use hnear combinations of Lorentz boosts to make 
the X- and F-components of both vectors proportional to each other with the same proportionality 
constant; then one rotates in the XY-plane to eliminate simultaneously the ^-components of both 
vectors. Alternatively, one can show that there are always suitable SL{2, M) transformations that 
simultaneously diagonalize two arbitrary matrices of the form (2.3). 
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Another consequence of the Einstein equations is that the functions /i,2(p) are not 
independent from each other: Both of them are solutions of the second order ODE 

/l,2 + 0Vl,2 = 0. 

Our bootstrapping algorithm to construct all solutions works as follows. With 
no loss of generality choose (2.9) and pick some arbitrary function /i(p). Then the 
first integral (2.6) dictates that /2 be a solution of 

/2/1 - /1/2 = 2j = const. (2.11) 

The constant j determines the angular momentum in (2.6) as J = (0,0,— j). The 
first order ODE (2.11) can be interpreted as Wronskian. Consequently, two linearly 
dependent functions /i oc /2 always require vanishing j, while non- vanishing j implies 
their linear independence. 

The next step determines the scalar field as a function of p by virtue of the 
Einstein equation (2.5a). (For later convenience we fix a sign ambiguity by assuming 
that (p is negative.) 

rr-H (2.12) 

Note that a reality condition on the scalar field implies a convexity /concavity 
property for the function /i. To be precise, /i has to be concave for positive values 
of /i and convex for negative values of /i. Thus, not any choice of the function 
/i leads to a suitable solution. Rather, we have to demand that in the range of 
definition of p the following inequality is valid. 

^ < (2.13) 

Ji 

This implies that /i cannot have any divergence at finite p. The same is true for /2. 

The next-to-last step is purely algebraic: The field equation (2.7) determines the 
potential ^ as a function of p. 

V = 2-i^(AA) (2.14) 

At this stage we know (p and V as functions of p. Since (f){p), when it is non-constant, 
is strictly monotonous within the range of definition we can, at least in principle, 
uniquely determine p{(f)). Thus, implicitly we know also the potential V{(f)). 

With the algorithm so far we can obtain all local solutions by scanning through 
all input functions fi compatible with the inequality (2.13). The last step is to relate 
the coordinates and used in the local construction above to time t and angular 
coordinate ip defined by their periodicity properties {t, (p) ~ (t, + 2n). They must 
be linear combinations of the coordinates (x^, x^), so that 



ax 



^ + bx^ ip = cx^ + dx^ (2.15) 
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Clearly, the above transformations generate solutions of the EOM with the same 
scalar potential. General solutions will then have the form 

= -f2{p){dxr + Mp){dx'r + TT^TT (2.16) 

/1(P)/2(P) 

identifying the angular and time coordinates according to (2.15). For static space- 
times we fix a = (i = 1 and 6 = c = so that (p = and t = x^ 

ds' = -Up) dt' + Mp) d^' + /^f (2.17) 

/l(p)/2(p) 

Some useful geometric quantities are the Ricci scalar 

R = -hh - /2A - \ hh = 2/1/2 0' - ^ /1/2 (2.18) 

the square of the tracefree Ricci tensor = R^,^ — | g^uR 

= \fifU' (2.19) 

and the Cotton tensor (see e.g. [5] for a definition) 

C,u = \hf2^'{5l5l + 5y^. (2.20) 

All solutions with j = or constant scalar field, = 0, are conformally fiat. 

We finally note that a given solution determined by the functions /i,2 permits 
to construct new solutions by taking linear combinations 

/i = A/1 + 5/2 /2 = C/i + D/2 (2.21) 

where A, B, C, D are some real numbers. Feeding the functions /i,2 into our algorithm 
then leads to new solutions analog to (2.16), since (2.11) is invariant under (2.21), 
with J = j {AD — BC). All these solutions have the same scalar field (2.12), but 
in general not the same potential (2.14). Thus, these solutions are in general not 
solutions of the same theory. A notable exception is A = D = 0, -B = C = 1, which 
we shall use in section 2.3 to construct solitons from black hole solutions. 

From the solutions (2.16), (2.17) it is clear that the function /i^2 should be 
chosen conveniently, depending on the application. "Conveniently chosen" refers 
to the sub-sector of solutions we are interested in — for instance, we may wish to 
restrict to asymptotic AdS solutions or to asymptotically fiat solutions, to solutions 
with or without horizons, etc. Moreover, certain conditions have to be imposed on 
the functions /i,2 to guarantee the absence of naked singularities. Imposing all these 
conditions from the very beginning of the algorithm is an important part of it, as 
this restricts the possible choices for the function /i and the integration constants 
appearing in the function /2. 
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2.2 General features — asymptotics, singularities, horizons, and centers 

In this section we analyze how to engineer solutions of EDG with the desired prop- 
erties. We assume throughout in this paper that spacetime is non-compact and that 
the asymptotic region is reached in the limit p — oo. This means that we pre-select 
in all our constructions the functions /i,2 such that this property holds. 

If we move from the asymptotic region towards the interior by considering smaller 
values of the "radial" coordinate p then eventually we may hit a singularity (curvature 
singularity, conical singularity, region with closed time-like curves or more exotic 
singularities); alternatively, we might not encounter any singularity but continue 
until we reach another asymptotic region in the limit p — )■ — oo. For most applications 
we are not interested in either of these possibilities. Instead, we would like spacetime 
to either contain a Killing horizon, shielding possible singularities (or superfluous 
asymptotic regions), or to contain a center, by which we mean a region around which 
spacetime looks locally like flat spacetime in polar coordinates close to the origin. 
In the former case an outside observer has access to the exterior region outside the 
Killing horizon, p > Ph, while in the latter case the whole spacetime is restricted to 
the semi-inflnite interval p G [po, oo), where p = po is the locus of the center. It is 
therefore of interest to establish a necessary condition on the functions /i,2 for the 
emergence of a Killing horizon or a center. 

Restrictions from regularity We consider the restrictions that smoothness around 
a chosen value of the radial coordinate imposes on solutions of the EOM. Assuming 
that the leading behavior at p = p is power law, two linearly independent functions 
satisfying (2.11) are 

/oc(p-p)" /cx(p-p)i-" (2.22) 

where a > 1/2 and we have ignored subleading terms at p ~ p. [For a = 1/2, 
/oc/log(p — p).] From (2.12), a pair /i, /2 constructed from linear combinations of 
the above two functions gives rise to the following behavior for the scalar fleld 

• y/a{l - a) 

= , ^ 

\P- P\ 

Reality of the scalar around p requires a < 1. Besides, unless a = l the scalar fleld 
develops a logarithmic singularity and, as can be seen from (2.18)-(2.20), there is a 
curvature singularity at the same locus. Having a equal unity is however not enough 
for smoothness at p. According to (2.12) we also need to demand 

/(p) = 0. (2.24) 

If this condition is not satisfled, including the flrst subleading term at p we have 

/ « (P - P) + ^/(P)(P - P)' / « 1 + /(P)(P - P) log(p - P) (2.25) 
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(2.23) 



which induces a curvature singularity. 

A similar analysis can be done in the asymptotic region. Regularity of the scalar 
field at infinity implies that /i, /2 are linear combinations of functions whose leading 
behavior at large values of the radial coordinate is / oc p and / oc 1. It leads to a 
constant value of the Ricci scalar at infinity, R — )■ — 1/1/2, and a vanishing Cotton 
tensor. Moreover, the metric asymptotes locally to (A)dS when fi and /2 grow 
linearly and the relative sign of both functions is (equal) opposite. If either /i or /2 
tends to a constant, the curvature vanishes at large p and the spacetime is locally 
asymptotically flat. In the following we shall restrict to these cases; namely we 
require the scalar to be bounded at infinity. An unbounded behavior of the scalar at 
infinity can be of interest for holographic applications, but we shall not address that 
situation in the present paper. 

Zeroes of the functions /i 2 Without loss of generality we choose /i to be asymp- 
totically positive. Taking fi as input in our algorithm, we obtain /2 integrating (2.11) 

/2 = x/i - 2j7i j Y (2.26) 
p 

where x ^md j are integration constants of the EOM, j having the interpretation 
of particle angular momentum in the reduced model (2.4). The second term on the 
right hand side of (2.26) provides a second solution, linearly independent from /i, to 
(2.11). If fi grows linearly at large values of the radial coordinate, this second term 
tends to a constant and vice versa. 

The inequality (2.13) implies that fi has a zero at a finite value of the radial 
coordinate, pi, if we assume regularity of /i for pi <p<oo at least up to its second 
derivative. We want to determine under which conditions /2 can have its first zero 
at p2>Pi- This is equivalent to search for the vanishing of 

X-2,/'| (2^27) 

P 

where we have introduced the upper integration limit R for concreteness; its only 
effect is to shift the value of the constant x- The behavior of the integral at pi 
depends on how fast /i approaches zero there. Let us consider first that it vanishes 
sufficiently fast to induce a divergence. The integral is a monotonously decreasing 
function. When /i oc p at large values of the radial coordinate, it is convenient to 
set R = oo such that the integral vanishes asymptotically. Hence (2.27) will have one 
and only one zero for any j/x > 0. By construction p2 is a single zero of /2 with 
vanishing second derivative. 

The sign of x determines whether the spacetime approaches AdS or dS asymp- 
totically. When x = the spacetime is asymptotically fiat, but there is no value of 



-9 - 



the radial coordinate where (2.27) becomes zero. An asymptotically fiat spacetime 
is also obtained for fi approaching a constant at infinity. In this case the integral 
in (2.27) diverges at large p to minus infinity, and that expression has one and only 
one zero for any j and x- 

When instead the function /i does not vanish fast enough at pi to generate a 
divergent contribution to the integral, there will be a maximum value of j/x for which 
(2.27) can have a zero at p2>Pi- At the limiting value the zero occurs precisely at 
pi. In the previous case, as j/x^^~^ we have p2— ^Pi- 
Let us summarize our results. Assuming the regularity of /i up to its second 
derivative and pi < p2, the function /2 has a regular zero at the point p2: a simple zero 
with vanishing second derivative. If pi is also a regular zero of /i the solution can 
be prolonged beyond this point, otherwise the lowest value of the radial coordinate 
where the solution is defined is po = pi- Let us stress the importance of the regularity 
condition on /i. If /2 has a simple zero with non vanishing second derivative, the 
function fi is finite but not differentiable at that point, see (2.25). As we have 
discussed this induces a curvature singularity, implying po = P2- Similar conclusions 
hold when P2 < Pi- Regular zeroes of the functions /i and /2 must alternate. [This 
as can be seen in the very simple example /i = cos p, /2 = sin p. Note, however, that 
this example leads to a linearly diverging scalar field at infinity.] When pi = p2 and 
j't^O, necessarily there is a curvature singularity at that locus and then po = pi,2- 

ADM gauge fixing Let us write the line-element (2.16) explicitly in terms of the 
time and angular coordinates 

ds^ = -AT^ dt^ + ^ + R^dip + dt) ' (2.28) 
hj2 

where 

R^ = a^f,-b^f, N^ = ^ N,= ^-^^^^ (2.29) 

Regions with negative contain closed time-like curves. Hence we shall always 
assume R^ to be positive in the asymptotic region. Defining a new radial variable 
as dr^ = dp^/R"^, the usual ADM form of the metric is recovered, with being the 
lapse function and contributing to the shift vector. 

The constants a,b,c,d were introduced in (2.15). The coordinate change (2.15) 
with the restriction ad—hc = \ provides a realization of the S'L(2,M) group of trans- 
formation of the matrix A^^, which locally preserves the line-element (2.2). A global 
rescaling of a, 6, c, d implements the multiplication of by a constant that, together 
with a compensating rescaling of the radial coordinate, also preserves Clement's pa- 
rameterization of the line-element. It will be convenient however to implement global 
rescaling directly on X, while the SL{2, M) group is realized in terms of the map be- 
tween local and physical variables (2.15). This choice has been already implemented 
in (2.28), where a,b,c,d are taken to satisfy the SL(2) condition. In order to avoid 
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redundancies in the local description, we need then to select a representative for each 
SL{2, M) orbit of the vector X. This is partly achieved with the gauge choice Y = 0, 
but it still allows for the family of solutions (7/1, /2/7) with 7GM. We will fix this 
remaining S'L(2,M) freedom as follows. When both /i^2 grow linearly at infinity, 
we require that /i — i- ±/2 asymptotically. When /i or /2 tend asymptotically to a 
constant, we will set this constant to ±1. 

Horizons A sufficient condition for the existence of a Killing horizon at a point 
ph of finite curvature is the vanishing of the lapse function N at that point. From 
(2.29) and the discussion above, this implies that /i, /2 or both must have a regular 
zero at p^. 

When j 7^ 0, the functions /i,2 are linearly independent and hence only one of 
them will vanish at ph- Let us assume that > at ph, so that we are not in 
a region of closed time-like curves. Reparametrizing the radial direction a.s y = 
a/2(p — Ph)/\j\ close to the zero, the line-element simplifies to 

ds^ = -e ^ dt^ + edy'' + Rlidip + N^hdtY + . . . (2.30) 

where the subindex '/i' indicates that the corresponding function is evaluated at ph- 
We have defined e = (— 1)* sign(j) with i = l{2) for vanishing /i(/2). When e = l the 
first two terms in (2.30) yield 2-dimensional Rindler spacetime, which is known to 
be the correct near horizon approximation of any non-extremal black hole Killing 
horizon (see e.g. [6]). When e = — 1 we obtain instead an inner black hole horizon 
or cosmological horizon. A consequence of (2.20) is that all Killing horizons have 
locally vanishing Cotton tensor. 

When j = 0, condition (2.11) implies the proportionality of /i 2. Therefore in this 
case both /i 2 vanish at a Killing horizon. Using then the coordinate transformation 

y = \l f{p- ph) we obtain 

ds' = -y'dt' + y' {dv + N,dt)' + j-^^ + ... (2.31) 

where according to our previous gauge choice /i,2 = /• The quantity A^^ is constant 
for vanishing j. For N^p = the metric (2.31) describes the Poincare patch horizon 
of AdS. 

Extremal black hole horizons are somewhat delicate with our gauge choices. A 
straightforward way to obtain them is from a 2-horizon configuration, performing 
suitable SL{2) transformations, and taking the limit where both horizons coincide 
in the end. In the near horizon approximation we have /i = 2(p — P4.) + 2£:p_|_, 
/2 = 2(p — p+), where e > is a small parameter that we send to zero eventually. 
We assumed here that at p = p+ > we have a black hole horizon [/2(p+) = 0] 
and at some slightly smaller value an inner horizon [/i(p+ — ^Ph-) = 0]. Our first 
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SL{2) transformation rescales the functions by (one over) e, /i = 2{p — p^)/e + 2p+, 
f2 = 2(p — p+)e. Our second SL{2) transformation is generated by (2.15) with 
a = c = 1, 6 = 1/e: — 1/2, d = l/e+1/2. Note that this transformation is singular 
at £ = 0, which is why so far we have to work at finite e. However, after these 
transformations we can smoothly take the limit e — ?• and obtain the near-horizon 
line-element 

2 dt^2(p-p+)^ dp^ , P+ aA"^ fo-io^ 

ds = + - r^ + 2p(d(^ dt) . 2.32 

P 4(p-p+)2 p 

Introducing r = a/2p leads to the extremal BTZ black hole in one of the standard 
coordinate systems, which is the near horizon limit of any extremal black hole in 
EDG. Note that close to p = p+ the first two terms in the line-element (2.32) form 
the expected AdS2 factor, see for instance [7]. 

Global restrictions We consider now to what extent the existence of a Killing 
horizon is compatible with different asymptotic behaviors. As before we assume fi 
to be asymptotically positive, so the different cases depend only on the asymptotic 
sign of the function /2, a = limp_>.oo sign (/2). As discussed after (2.25), a cannot be 
zero for asymptotically regular solutions, so we need to discuss only two cases. 

0" = + This covers asymptotically locally AdS and flat spacetimes. In order to have 

positive at the outermost horizon, the function /2 must vanish there, implying 
e = l. Behind the horizon we may encounter: i) a, curvature singularity where 
/i,2 vanish; ii) a. ring of curvature singularities of the type (2.25), where only /i 
vanishes while is positive; Hi) a regular zero of /i with e = — 1 producing an 
inner Killing horizon. Notice that an inner horizon is only present in the non- 
generic situation in), otherwise there is a curvature singularity. This agrees 
with the results of [8], obtained studying matter fluctuations around a BTZ 
black hole. For static spacetimes, the ring of singularities in case ii) collapses to 
a point, while the regular zero of /i in case Hi) gives rise to a conical singularity. 

a = — This covers asymptotically locally dS and asymptotically locally fiat cosmolog- 
ical solutions. The outermost Killing horizon necessarily satisfies e = — 1, and 
consistently has the interpretation of a cosmological horizon. Inside the cosmo- 
logical horizon it turns out not to be possible to find a second Killing horizon 
without reaching before a region with R^ negative and thus closed time-like 
curves. Hence this rules out the existence of regular black hole geometries in 
asymptotically dS spacetimes in the context of EDG. 

We will not have more to say in this paper about asymptotically dS or fiat cosmo- 
logical solutions (see [9]). Hence from now on we focus on asymptotically positive 
functions /i,2 as input for our algorithm. 
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Centers A necessary condition for the existence of a center at po is the vanishing 
of the radial function R while the lapse function remains finite. This can only 
be achieved if either a or 6 (but not both) vanish. Requiring /i^2 and to be 
asymptotically positive necessarily selects 6 = 0. Hence centers are associated with 
regular zeroes of /i. To reduce clutter we focus on static spacetimes. (Stationary 
spacetimes can be treated analogously.) Changing again the radial coordinate to 
y = a/2(p — Po)/|j|) the line-element close to such a zero simplifies to 

ds^ = -f2cdf + ^y^dif' + dy-' + ... (2.33) 
he 

with /2c = /2(po)- Smoothness requires furthermore the condition 

f2c = f. (2.34) 

If this condition is violated there is a conical defect at p = po. 

Notice that the relation (2.34) together with positivity of fic and the condition 
(2.11) imply that j must be negative for solutions without conical singularity. By 
contrast, static solutions with a horizon and a regular zero of /i discussed in the 
previous paragraph have positive j. Thus, there is always a conical singularity hidden 
behind a horizon in these cases. 

2.3 Constants of motion vs. parameters of the model 

The algorithm above involves a number of free constants and provides a potential V 
as an output that depends in general on these constants. However, once a potential 
V is constructed it is also of interest to obtain further solutions to the EOM with 
the same potential. These new solutions differ from the previous one by the values 
of certain constants of motion. 

Counting of constants of motion The Einstein equations (2.5a) and Klein- 
Gordon equations (2.5b) lead to eight integration constants. The Hamilton con- 
straint (2.5c) relates them, so that we have seven independent constants of motion. 
(Time derivatives of the Hamilton constraint do not lead to new relations between 
these constants of motion.) Three of these constants are contained in the 5*1/(2, M) 
transformations of the matrix A^,y in (2.3). Actually, two of the SL{2, M) transforma- 
tions have been used to eliminate the variable Y in the matrix (2.3). Since the action 
(2.4) does not depend explicitly on the radial coordinate, shifts of p generate new 
solutions. This freedom can be parameterized by the value of po. Another integration 
constant is the conserved particle angular momentum parameter j. The remaining 
two constants of motion, denoted by (p±, are coefficients of the two scalar modes and 
appear implicitly in the specification of the function /i, the starting point of our 
algorithm. We shall see explicitly the emergence of these constants when discussing 
asymptotically AdS solutions in sections 4 and 5. 
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Co-dimension 2 family of solutions Contrary to SL{2,'R) transformations and 
radial shifts, our algorithm does not instruct us how solutions must depend on j 
and (j)± in order for the scalar potential to be independent of them. Given /i we 
can construct /2 for any value of j from (2.26). But these solutions will not relate 
to the same potential unless /i carries an a priori involved dependence on j, which 
we do not know how to implement systematically. There is, however, a simple way 
of constructing a one parameter family of solutions to the same V{(f)) which moves 
in the interesting subspace of constants of motion {j,4>±}. It is generated by the 
rescaling of /i,2 together with a compensating rescaling of the radial coordinate 

fiM ^ ^fiM^) M ^ <PipM (2.35) 

Here /i and /2 are assumed to solve the EOM for some j and (f)±. Then 

j Vtj ^ fi^±/2(/.± (2.36) 

where A± are the scaling dimensions of (j)±. (Also po ^ ^po-) As already mentioned, 
transformation (2.35) is a symmetry of Clement's parametrization of the line-element, 
and hence of the particle model (2.4). This guarantees the scalar potential to be Q- 
independent. Alternatively, the fact that both the scalar field and the potential as 
function of the radial coordinate (2.14) depend on p and Q through p/Q, leads to the 
same conclusion. Motivated by the convenience to use Q as one of the constants of 
motion, we change parametrization in the space of solutions from {j, to {fl, ^, A}, 
where 

e = j/n A = (2.37) 

Analytic expansions As implied by (2.12), smooth solutions with a Killing hori- 
zon or a center require fi^2 to be at least twice differentiable down to its locus. 
In order to carry on a general analysis, we further assume now all functions to be 
analytic for p>po- Around a generic radial point p, we have 

/l=2^an /2 = 2^fen <P = }_^<Pn (2.38) 

n=0 n=0 n=0 

The conservation law (2.11) implies 

aoh - aibo = 2^ (2.39) 
- 02^0 = (2.40) 



which fixes the coefficients 6„>o in terms of ^, am<n and Bq. Similarly, the defining 
equation for the scalar field (2.12) determines the coefficients (pn>o in terms of am<n+2- 
The potential as a function of the radial coordinate is given by 

^ = 2 - i E - 1) E ^rnbn- J^'Z (2-41) 

n=2 m=0 
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On the other hand, assuming that the potential V{(j)) is also an analytic function of 
(j) around (^(p) =Co yields 

V = V(0o) + (0 - hWih) + ^^^^ V"(0o) + . . . (2.42) 
Equating the two expansions (2.41) and (2.42) we obtain to leading order 

V{4>o) = 2 - ^ (ao62 + ai6i + aa^o) (2.43) 

This relation determines the coefficient 02 in terms ao,bo,ai,(l)Q and Comparison 
to higher orders fixes a„>2 as functions of the same parameters. The parameters oq, 
bo and Q are not independent from each other; instead, one of them can be fixed in 
terms of the others. Consistently, we recover that the space of solutions of EDG in 
the gauge y = is five-dimensional. 

Static black holes and their associated solitons We focus now on the inter- 
esting case of having a black hole horizon at p. As discussed in the previous section, 
this requires to impose Bq = 0, while the regularity condition (2.24), viz. 62 = 0, is 
automatically satisfied by (2.40). Using the freedom to redefine Q we can set the 
product Oo^i to any desired value. For convenience we choose 

aobi = 8 . (2.44) 

Relation (2.39) fixes then ^ = 4, eliminating a parameter with respect to the generic 
case. With this choice the relation (2.43) and higher orders determine a„>o in terms 
of bi/ttQ and 0o- Allowing for radial shifts, we are left with four integration constants. 
Hence black hole solutions span a co-dimension 1 surface in the space of solutions to 
EDG with a chosen scalar potential. 

Every fi-family of /i,2 functions giving rise to static spacetimes with a (non- 
extremal) Killing horizon, has an associated smooth solitonic solution with a center. 
(Similar conclusions hold for stationary solutions.) This solution is obtained by 
exchanging the roles of the input functions 

fls = f2 /2s = /l (2.45) 

and demanding the absence of a conical defect. The absence of a conical defect, 
condition (2.34), specifies uniquely the value of Q to 

ns = ^ (2.46) 
lb 

where we have used (2.44). Therefore, for any (n + l)-parameter family of static 
black holes there is a smooth n-parameter family of static solitonic solutions with a 
center. The simplest example is the 1-parameter family of static BTZ black holes, 
which has a 0-parameter 'family of solitons', namely the global AdS solution. It is 
suggestive to consider these solitonic solutions as possible ground states of a given 
theory. To decide which geometry is the ground state we have to calculate the free 
energy, which we shall do in section 6. 
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Rotating black holes Locally there is no difference between rotating and static 
blaclc holes in three dimensions. Indeed, all BTZ black holes [10] — rotating or static 
— are locally AdS and thus locally equal to each other. Rotating black holes differ 
from static ones by their periodicity properties of the coordinates x^. For the 
same functions /i and /2 that produce static black holes with the choices a = d = 
1, 6 = c = in (2.15), we can pick other values for a, b, c, d to obtain rotating 
black holes. In parlance of canonical analysis the SL(2) transformations that we 
exploited for a local gauge fixing may fail to be gauge transformations globally, since 
their canonical boundary charges can be non-zero. When studying global solutions 
we should therefore undo the gauge fixing that led to the Ansatz (2.10). Of the 
three independent SL{2) transformations, one leads to a rescaling (7/1, fill)-, while 
another one leads to a rotation of the asymptotic metric. We are not interested here 
in either of these transformations, but instead want to consider only those SL{2) 
transformations that keep fixed the asymptotic metric. These transformations are 
generated by a single parameter rj. 

For black holes in AdS, with our gauge choice limp_^oo /1//2 = 1? the transfor- 
mations are 5*0(1, 1). 

t = cosh 77 — x^ sinh?7 (p = —x^ sinh?7 + x^ cosh?] (2.47) 

Static solutions (2.17) are then the 77 — )■ limit of a larger set of stationary solutions 
(2.28) with the identification (2.47). Similarly, for asymptotically fiat solutions with, 
say, limp^oo /i = 1? linip-5.oo /2/p = 2, the relevant transformations are parabolic, 
t = x^, (f = —x^ 7] + x^.^ Finally, for dS the transformations are 5*0(2); but recall 
that in EDG there are no black holes in dS, as we have shown in the previous section. 

It makes sense to consider the parameter 77 as one of the constants of motions, 
since it parametrizes the black hole angular momentum. Thus, globally the total 
number of potentially relevant constants of motions is four: {77, fi, ^, A}. 

Concluding remarks We have now a clear understanding of the role of nearly all 
the seven constants of motion: one is fixed trivially by choosing the value of the radial 
coordinate at the center /singularity, po- Two, related to SL(2) transformations, are 
fixed by choosing the boundary metric. The remaining four are physically more 
relevant, with the following meanings. The parameter 7/ (related to residual SL{2) 
transformations) generates rotating solutions from static ones. The parameter Q 
generates a 1-parameter family of solutions for the same scalar potential, in general 
with different masses. The parameter A will be discussed in more detail in section 
4. We finish this discussion speculating on what type of deformation can move away 

^Notice that allowing 6 7^ in this case will induce negative at infinity and thus closed 
time-like curves. This provides an alternative argument for disregarding one of the three SL(2) 
parameters in this type of asymptotically flat solutions. 
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from geometries with a Killing horizon while keeping the same scalar potential, which 
helps us to identify the role played by the seventh parameter ^. 

In the past subsection we have explored the construction of asymptotically flat 
or (A)dS solutions, taking as input a function /i positive at large values of the radial 
coordinate and at least second order differentiable up to its first zero. The inequality 
(2.13) forces /i to vanish at some finite radial point po- A sufficient condition for the 
existence of a Killing horizon, which would hide the generically singular behavior at 
Po, is the divergence of the integral fdp'f^'^ as p tends to po. Therefore, under the 
mentioned differentiability assumption of /i, a necessary criterion for the absence of 
a Killing horizon is the finiteness of the previous integral at po- 

The no-Killing-horizon criterion is fulfilled when fi = ai_Qr2"(p — po)^~" with 
1 > a > 1/2 close to po. Varying a can be naturally interpreted as a deformation from 
the case with a horizon if at the same time /2 = &afi^~"(p— Po)"- A horizon at po can 
be obtained when a equals unity, while moving to smaller values of this parameter 
leads to geometries with a naked singularity.^ Moreover, for the proposed /i and /2 
the relation between ^ and a is 

^ = 8a - 4 (2.48) 

when, analogously to (2.44), we fix ai-Jja = 8 by an appropriate choice of Vt. The 
complete form of the functions /i^2 should be determined by keeping the scalar poten- 
tial unmodified. Since in the solutions we are considering cj) diverges at po [see (2.23)], 
this question concerns the form of the potential at large values of the scalar field. 
The discussion above suggests that varying ^ in general leads to singular solutions. 
We analyze next a very simple example where this picture is realized. 

2.4 Solutions with vanishing potential — gravity plus kinetic energy 

For V = all solutions can be constructed in closed form (see for instance [11]). This 
is one of the rare cases where our algorithm can be inverted and taking the potential 
as input explicitly leads to the functions /i,2 as output. Integrating (2.7) with the 
left-hand-side set to zero obtains 

= 4(p2 + cip + C2) (2.49) 

where Ci,2 are integration constants. Depending on the sign of the discriminant 
A = cf — 4c2 we obtain different cases. We discuss explicitly the case of positive 
discriminant. 

With no loss of generality we fix ci = by exploiting shifts and define C2 = — 
Integrating the Klein-Gordon equation (2.5b) twice yields 

^Sincc /2 vanishes at po faster than /i, expression (2.27) must become zero at po- Given that 
(2.27) can only have one zero for p>po, the associated singularity will be naked. 



-17- 



where (f)± are integration constants. The value of 0+ is bounded by < 1. 

Violating this bound leads to a change of the sign of the discriminant. Integrating the 
Einstein equation (2.5a) and imposing the Hamilton constraint (2.5c) then establishes 



n,^2V^^(^-±^f^ (2.51) 



which introduces one more integration constant. 

The Ricci scalar (2.18) and the Cotton tensor (2.20) are given by 

R- 6+ ''^'^ C -^A^lEMl (252) 

with other components of C^i, vanishing (besides the obvious Cio = Cqi). These so- 
lutions are always asymptotically AdS. If 0^ = or 0^ = r2^ they are also conformally 
flat. At = there is a curvature singularity for non-zero 

When 0+ vanishes we obtain vacuum solutions and the function (2.51) simplifies 

to 

/i = 2(p + fi) /2 = 2(p-fi) (2.53) 

with a constant scalar = 0_. The result (2.53) explains why we chose to fix ao^i = 8 
in section 2.3, since we obtain fi\p^n = 4f2 and f2\p^n = 2(p — fi), so oq = 4 and 
bi = 2. This family of solutions should reproduce the BTZ geometries [10]. Indeed, 
redefining the radial coordinate as r = ■\/2(p -|- f2) and identifying coordinates as in 
(2.47) the ADM form (2.28) of a BTZ black hole is recovered, with mass and angular 
momentum 

M = 87r cosh 2r/ J = 87r 1] sinh 2r/ . (2.54) 

Extremal BTZ emerges as a double scaling limit 1] — 0, |?7| — oo. An alternative 
construction of extremal BTZ is possible along the lines of section 2.2 and leads to 
the line-element (2.32). 

When is non- vanishing, the behavior around p = is of the form 

/i =2^+"fi"(p-fi)i-° /a = 22-"fii"" (p - fi)° (2.55) 



with a = (l + yl — 0^/r2^)/2; thus l>a>l/2. We then recover precisely the result 
(2.48). These solutions realize the scheme discussed above. Solutions with a Killing 
horizon span a co-dimension one subspace in the space of constants of motion. The 
parameter ^ moves away from geometries with a Killing horizon into others with 
naked singularities. Moreover, the value of ^ is only determined by the exponent a 
characterizing the singularity. 

The cases of negative and vanishing discriminant lead to complex solutions. 
There are no non-trivial real solutions for the metric and the scalar field in these 
cases. In conclusion, all non-trivial ^ = solutions with real metric and real scalar 
field have naked singularities, in accordance with the "no-hair" property of black 
holes. The only real and regular solutions for vanishing (or constant) potential are 
vacuum solutions. 
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3. Asymptotically flat solutions 



— black holes with scalar hair 



In this section we present an explicit example of locally asymptotically flat black 
holes with scalar hair. We define locally asymptotically flat solutions by the property 
that Ricci scalar and Cotton tensor asymptote to zero for p — )■ oo, and the metric 
asymptotes to the fiat metric in the same limit. The scalar field and its p-derivative 
then automatically are finite in the limit of large p [see (2.12)], and the potential 
goes to a constant [see (2.7)]. 

Before starting let us briefly summarize a key result from section 2, which leads 
to the following no-go theorem: there are no globally asymptotically fiat hairy black 
hole solutions in EDG. Indeed, in order to obtain an asymptotically fiat spacetime 
we need that either /i or /2 tend to a constant at infinity. We saw in 2.2 that the 
existence of a Killing horizon is only possible when the bounded function at infinity 
is /i. This implies that the space at infinity contains a circle of finite radius, and 
thus can only be locally asymptotically flat. 

Let us therefore consider locally asymptotically flat black holes with scalar hair. 
A simple choice for a function /i approaching a positive constant asymptotically is 

/i = l--. (3.1) 
P 

Reality of the scalar field requires r2>0. From (2.26) we have 

/-"V.(-4-2«[;r^-^2.og(i^-l)]) (3.2) 

The choices above with the identification x° = t, = lead to a line-element that is 
asymptotically Rindler2 x S^. The function /i has a simple zero with non- vanishing 
second derivative. This leads to a curvature singularity at its location p = fl 

i? = -6eiog(^-l) + ... (3.3) 

The singularity is always covered by a Killing horizon, since the term in parenthesis 
in (3.2) has one and only one zero for any A, Q and ^. The scalar field is obtained 
integrating (2.12). 

= 00 + V2tt - 2V2 arctan - 1 (3.4) 

It has a divergent kinetic energy at the singularity, but it is regular otherwise. Thus, 
these solutions are explicit examples of asymptotically flat black hole solutions with 
non-trivial scalar hair. 

The associated scalar potential is derived as a function of the radial coordinate 
from (2.14). The inversion p = p(0) is easy in this case, and we obtain 

V = 2 + ^ sin^0+ (y4-2^ + 4^1ncot2 0) sin^0-| {A + i^ + i^lncot^ 4)) sin^ (3.5) 
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V 

2.5 r 




Figure 1: Potential for the exact solution (3.5) plotted as function of sini?!> for the choices 
Q = 1, A = 2, ^ = i and (/)o = 

where we have defined 0= (0— 0o)/"\/8- As expected the potential is independent of 
the parameter Q. It does, however, depend on A and C,- In subsection 2.3 we have 
considered C, as one of the constants of motion that parametrize solutions of EDG 
with the same potential. Preserving V required a very non-trivial dependence of 
both /i 2 on ^, codified in (2.39), (2.40) and (2.43). Since our algorithm does not tell 
us how to solve those equations, we have constructed the present example based on a 
^-independent function /i. When ^ is implemented as constant of motion, we argued 
that it generically moves away from solutions with a Killing horizon into ones with 
naked singularities. This is avoided in the present case by the explicit dependence of 
the potential on this parameter. 

The potential (3.5) is plotted in figure 1. In addition to A and ^, it depends on 
the integration constant of our algorithm 0o through the redefined field </>. It has a 
local minimum at (f) = (j)Q, two global minima lying symmetrically around 4' = 4>o and 
diverges if sin0— )-±l. It is important to recall that only for a constant potential the 
quantity 0o is a constant of motion in the sense considered above. 

To the best of our knowledge solutions (3.1)-(3.5) are the first explicit examples 
of locally asymptotically flat black hole solutions with non-trivial scalar hair in 3- 
dimensional EDG. The evasion of no-hair theorems has been possible, because these 
theorems assume typically V > 2 (see [12] and Refs. therein), which is easy to violate. 
As evident from figure 1 the potential (3.5) violates the inequality V > 2 in some 
range of the coordinate p. Physically the evasion of no-hair arguments is possible 
due to the presence of a minimum below the asymptotic value of the potential, 
limp_>.oo V = 2; see e.g. the discussion in [13]. 
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4. Asymptotically AdS solutions and scalar asymptotics 

From the line-element (2.16) we see that in order to obtain an asymptotically AdS 
metric with a Minkowski boundary metric the functions /i,2 must tend to the same 
linear function at large p, namely 

/i,2 = 2p+--- (4.1) 

where the dots stand for subleading terms and the overall factor is fixed by requiring 
unit AdS radius. This fixes the boundary metric and so it corresponds to a Dirich- 
let boundary condition for the metric. However, in order to properly define the 
variational problem for EDG we need, in addition, to prescribe suitable boundary 
conditions for the scalar field </>. In the current section we discuss a simple class of 
boundary conditions for the metric and scalar fields that preserve the AdS asymp- 
totics and which suitably generalize the Brown-Henneaux boundary conditions for 
pure gravity [14] . A complete analysis of the most general boundary conditions for 
EDG gravity that lead to AdS asymptotics is performed in section 5. 

4.1 Boundary conditions — Brown— Henneaux and beyond 

The Brown-Henneaux boundary conditions for vacuum 3D Einstein gravity with a 
negative cosmological constant correspond to the Dirichlet boundary condition (4.1) 
for the metric. However, this leading asymptotic behavior in fact determines through 
the EOM certain subleading terms as well. In particular, for vacuum Einstein gravity, 
the Brown-Henneaux boundary conditions are'^ 

/i,2 = 2p + Ci,2 + o(l) (4.2) 

where Ci^2 are arbitrary constants, whose difference determines the parameter j 
(2.11). Introducing a canonical radial coordinate through 

J V flip) Up) V 8p J 

where c 7^ is an arbitrary integration constant, the line-element (2.16) becomes 

= ^ + [glf + g'ff + g-^ dx^ dx^ (4.4) 

where 

= 2c-2diag(-l, 1) g('^ = diag(l, 1) g'^^' = o(l) (4.5) 

^The notation F{x) = o{f(x)) as x ^ xq means limj;^.^;;, jr^ = 0. 
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The expansion (4.4) contains precisely the leading background term and the sub- 
leading state-dependent term g^"^^ that appear in the original Brown-Henneaux anal- 
ysis [14]. However, neither g^'^\ nor g^'^^^ are part of the specification of the bound- 
ary conditions. In the case of vacuum Einstein equations one in fact finds that 
^rest _ Q^j^-^^j^ Iq^^i ^j-^is will uot in general be the case for EDG. 

For EDG instead, Brown-Henneaux boundary conditions correspond to the 
asymptotic expansions 

/i,2 = 2p + Ci,2 + 0(p-^) (4.6) 

for some unspecified e > 0. Inserting these into the equation (2.12) for the scalar 
field leads to the corresponding asymptotic form of the scalar. Namely, 

0(p) = O(p-^) (4.7) 

This asymptotic behavior in turn determines the mass term in the expansion of the 
potential around = 

y=lm2 02 + --- (4.8) 

where the scalar mass is determined in terms of e as 

m^ = ~{l-e^) (4.9) 

It follows that EDG can admit asymptotically AdS solutions provided the parameter 
e introduced in (4.6) lies in the range 

< e < 1 (4.10) 

where the lower limit corresponds to the BP bound = — 1 and the upper one 
to = 0. However, once the scalar mass is determined so are the two possible 
asymptotic behaviors of the scalar field, namely 



A 



± . 



0(p) ~ 0{p~-) (4.11) 
with 

A± = l±£. (4.12) 

A+ then corresponds to the Brown-Henneaux boundary conditions (4.7), while A_ 
corresponds to a more general class of boundary conditions. 

Since the scalar field is determined from the input function /i(p) via the ODE 
(2.12), in order to generalize the boundary conditions and capture both possible 
asymptotic behaviors of the scalar we need to suitably generalize the asymptotic 
expansions (4.6). This can be done systematically and in complete generality. We 
do this in section 5. For now, we discuss the simplest possible generalization of (4.6), 
which takes the form 

/i = 2p + An^-' + Bn logp + c^] + ■■• (4.13) 
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where A, B, C and 7^ are constants, and the dots stand for terms suppressed at 
large p. Moreover, we assume < e < ^ for now. 

Given the asymptotic expansion for the input function /i(p), the constraint 
(2.11) can be integrated to obtain 

/2 = /i-e^^ + --- (4.14) 

with j = ^Vt. Moreover, inserting the expansion (4.13) into (2.12) we get 

which can be integrated to obtain 

= 0„p^^-/2 + 0+p^^+/2 + --- (4.16) 

where 

ct^. = an^-" ct>^ = ^^n^-" (4.17) 

aA_A_|_ 

and we have defined 

/ 2Ae 

a^^—. (4.18) 

We see that the input function (4.13) leads to the correct asymptotic form of the 
scalar with both possible asymptotic modes. It is therefore a valid generalization of 
the Brown-Henneaux boundary conditions for < £ < 1/2. 

The parameter B is proportional to the integration constant A defined in (2.37) 

B = A_A+«^A (4.19) 

The proportionality constant is set by the exponent e, which fixes the scalar mass 
and thus is a characteristic of the potential. We have argued in subsection 2.3 that 
starting with a regular solution and varying ^ while keeping fi, A and the scalar 
potential fixed, in general leads to singular solutions. Since mapping to the dual 
CFT is only established for regular solutions, we will not consider the dependence of 
solutions on ^. On the other hand, implementing B as integration constant would 
require a very non-trivial dependence on this parameter of the /i piece unspecified in 
(4.13), on which we will not enter. As extensively discussed only the (^-dependence 
of solutions is immediate to implement, which we have done in (4.13). 

The usual AdS/CFT dictionary interprets 0_ as a source term in the dual CFT. 
According to this, the 1-parameter family of solutions generated by Vl would move 
among boundary field theories with different couplings. However, since both modes 
of a scalar with mass in the window (4.9), (4.10) are normalizable, we can choose 
not only 0_ as a source (Dirichlet boundary conditions), but also 0+ (Neumann 
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boundary condition) or any combination of the modes and 0+ [15]. Especially 
interesting boundary conditions are those that allow to vary Q without modifying 
the boundary field theory. Clearly this is achieved if we impose 

0+ - A^^+Z"^- = (4.20) 

with A fixed. This boundary condition corresponds [15] to subject the CFT associated 
to Neumann boundary conditions to a marginal multitrace deformation 

2 

AL = \O^Z (4-21) 

where 0a_ is the single trace operator dual to 0+. Black hole families generated by 
Q provide then the holographic dual for certain thermal states in such a deformed 
CFT. In section 6 we will discuss their thermodynamic stability, as well as the role 
of the associated soliton geometry (2.45), (2.46) [16,17]. 

Since we want to study the thermodynamics of the dual field theories, we need 
to address the holographic renormalization of our backgrounds. Although practical 
reasons highlight marginal boundary conditions, any other ones are a priori of equal 
relevance. Notice that the stability of a solution depends crucially on the chosen 
boundary conditions. We will start below the analysis of holographic renormal- 
ization considering Dirichlet boundary conditions. Neumann and Mixed boundary 
conditions will be treated in section 4.3. 



4.2 Holographic renormalization 

The bulk action for field configurations solving the EOM can be rendered finite by 
the addition of appropriate boundary counterterms 



R 



{d^f - V{cf>) 



d^XA/— 7 



K - U{<j)) (4.22) 



where the subscript '+' indicates that we choose Dirichlet boundary conditions for 
the scalar. K is the trace of the extrinsic curvature on the boundary 9, defined as 
a. p = const, surface for asymptotically large values of the radial coordinate, and 7 
is the induced metric on d. For solutions corresponding to an input function of the 
form (4.13) the counterterm function U{(f)) is given by 



W) = l + — 



(4.23) 



We will derive this result in section 5, where the form of U{(j)) for the most general 
asymptotically AdS solutions is also derived. However, we emphasize the results we 
give in the following for the renormalized action and one-point functions are in fact 
valid for the most general asymptotically AdS solutions that depend only on the 
radial coordinate, as can be explicitly confirmed using the counteterms derived in 
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section 5. Moreover, it should be stressed that the counterterms in (4.22) are not 
sufficient for regularizing general solutions, but only those which solely depend on the 
radial coordinate. Since in this paper we are considering stationary, axisymmetric 
solutions of EDG, it is enough for our purposes. 

Taking into account all boundary terms, the on-shell action is 



Vol 



I EOM 



5p(/i/2)|p, + lim 



5p(/l/2)-4t/(0)v^ 



Vol(|j|+A_(A+-A_ 



(4.24) 
(4.25) 



with Vol = Jdx^dx^. The quantity ph denotes the locus of a Killing horizon or 
center, which is the only natural lower bound of the radial integration domain. The 
discussion in section 2.3 shows that when the asymptotic signs of /i and /2 coincide 
the reduced model angular momentum j is positive for metrics with a Killing horizon 
while negative for those with a center. Hence the ffist term on the right hand side 
of (4.24) defines the absolute value of j on solutions of the form (4.13). 

The variation of the full on-shell action (4.22) with respect to the induced metric 



ST, 



I EOM 



-7 



(ir-f/(0))7i, Sf^ 



(4.26) 



determines the holographically renormalized Brown- York stress tensor, i.e. the renor- 
malized one-point function of the stress tensor of the dual field theory, namely 



5r. 



gf0)=2r, 



(4.27) 



where we have set the constant c = 1 in (4.5) for solutions of the form (4.13) so that 
lij^dfj* P + • • • ) with glj^ = 2'r]ij. Evaluating this we obtain 

(T,,)+ = J 5,j + A_(A+-A_)0+0_ r/,, (4.28) 



whose trace is 



9 



(0)ii 



(T,,)+ = A4A+-A. 



(4.29) 



Notice that if either vanish then the stress tensor is traceless. If additionally 
the functions fi are proportional to each other then j vanishes, the spacetime is 
conformally flat and the energy momentum tensor vanishes completely. 

Another interesting object is the response with respect to variations of the scalar 
field, which provides the expectation value of the operator dual to the 0_, 0+ com- 
bination that we chose as source. Dirichlet boundary conditions yield 



I EOM 



- / d^X. 
Id 



-1 



n^d^<p + 2f/'(0) 



5<p 



(4.30) 
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where is the outward pointing unit normal vector to the boundary d. The variation 
of the scalar field associated with the mode 0+ gives an asymptotically vanishing 
contribution to (4.30), and therefore only 50 = (50_ p-^-/2 is relevant. The vacuum 
expectation value of the field theory operator dual to 0_ then reads 

(0..) = -^^ = -(A.-A_)^, (4.31) 

The ratio between the vacuum expectation value of the trace of the stress tensor and 
the scalar response is a 0+-independent constant 

= -A_0„ (OaJ (4.32) 

This is the usual trace Ward identity. 

4.3 Beyond Dirichlet boundary conditions 

We will analyze now Neumann and Mixed boundary conditions. We denote as 
i7(</)_, 0+) the combination of the two scalar modes that one decides to take as 
source. We can have 



Dirichlet J 
Neumann J' 
Mixed J 



(f)- (4.33a) 
(A+-A_)(/.+ (4.33b) 
(A+-A_)0+ - >V'(0_) (4.33c) 



where VV(0_) is an arbitrary function. The field theory operator sourced by J has 
conformal dimension A_|_ in the case of Dirichlet boundary conditions and A_ in the 
cases of Neumann or Mixed boundary conditions. 

In order to have a well posed variational principle, the change of the action 
under variations of the fields must be proportional to the variation of the sources 
(modulo the EOM). Indeed the sources are naturally to be held fixed for a dual 
field theory interpretation. This justifies the claim that F^, as defined in (4.22), 
corresponds to Dirichlet boundary conditions. Namely, its variation is proportional 
to (50_. Similarly, it follows that in order to impose Neumann or Mixed boundary 
conditions the following boundary term must be added to (4.22): 

Syv = - [ d2a;V-5(o)[(A+-A_)0+0_ + W(0_) - 0_>V'(0_)1 (4.34) 



The Neumann case is obtained by setting W = 0. Adding this extra boundary term 
we find that the change of the on-shell action under variations of the scalar field is 

5r_ = 5{T+ + Sw) = - j d''x^-g(^) (P^5J (4.35) 
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which is proportional to 5j7 as required. An immediate conclusion we draw from 
this identity is that the scalar one-point function for Neumann or Mixed boundary 
conditions is given by 

(C?A_) = -^i=^ = 0- (4.36) 

The one-point function of the stress tensor for Neumann or Mixed boundary 
conditions is given by [18] 

(T,,)_ = (T,,)+ - (W(0-) + 0- J) gf^ (4.37) 

The trace of the stress tensor yields 

{T\)_ = -A+0_ J + A_0_ W'(0_) - 2>V(0_) (4.38) 

Notice that for marginal boundary conditions (4.20), (4.21) we have 

W(0_) = A ^-^^;-^-^ 0^ (4.39) 

and the trace Ward identity simplifies to 

(n)_ = -A+J(OA_). (4.40) 

This is exactly of the form (4.32), except that expectation values and sources have 
been changed and the dual scalar operator in this case has dimension A_. 

We conclude this section with a couple of observations. Rewriting the boundary 
term (4.34) as 

Sw = - [ d^x^/^^LJ + W(0_)1 (4.41) 



and using the second equality in (4.36), it follows that the effective action for the 
scalar expectation value 0_, given by the Legendre transform, takes the form 

r[0_] = T_- Jd^xJ^J^=T+- J d^x^/^) >V(0_) . (4.42) 

Hence the function >V(</)_) corresponds to a multitrace deformation of the theory 
corresponding to pure Neumann boundary conditions. 

A second remark concerns the relation between boundary conditions and regu- 
larity of the solutions. States in the theory, i.e. regular solutions, corresponding to 
Mixed boundary conditions satisfy 

J = (A+-A_)(/)+->V'((/._) = 0. (4.43) 

This imposes a relation between the two modes 0_ and 0+. However, the condition 
that the solutions be also regular imposes another condition, which can be written 
in the form 

(A+ - A_)0+ - 9^_>Vreg(j, 0-) = 0, (4.44) 
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for an a priori different function VVreg(j, </>-)• Combining this with (4.31) and (4.42) 
we get 



up to an unphysical constant. This is the full effective potential for the VEV of the 
dual operator Oa_- Smooth EDG solutions correspond to extrema of this effective 
potential, i.e. they must satisfy (Wrcg — W) = 0. 

In section 2.3 we discussed the criterium for regularity of solutions in term of the 
constants of motion {Q, ^, A}. We showed that ^ is fixed by regularity. Moreover, 
defining Q to satisfy (2.44) lead to a constant value of ^ on regular solutions. This 
convenient choice was possible due to the freedom to redefine f2 by a multiplicative 
factor that can depend on the other constants of motion, which is manifest in (2.35). 
In this section we have performed an analysis of solutions of EDG based only on their 
asymptotic properties, while imposing (2.44) requires the knowledge of the interior 
geometry. Hence for the parameterization (4.13), the regularity condition will instead 
read 



This relation is equivalent to (4.44). 

Once the regularity condition is imposed, the smooth solutions are described by 
the parameters {Q, A}, or equivalently {j, 0_} (and globally in addition by the rota- 
tion parameter rj). If we impose marginal boundary conditions (4.39) the parameter 
A is fixed for all solutions compatible with these boundary conditions. For other 
boundary conditions A varies on the solutions. 

5. Asymptotically AdS solutions — general discussion 

In this section we generalize the discussion of the previous section by extending 
the range of allowed values of e in the asymptotic expansion (4.13) to the whole 
interval < e < 1. In particular, we determine the most general form of the input 
function /i(p) such that the theory admits AdS solutions and can be holographically 
renormalized by local counterterms. On the way we determine the general form of 
the counterterms in (4.22), as well as the most general asymptotic form of the scalar 
potential. 

5.1 Most general asymptotically AdS solutions 

It must be emphasized that we are not interested in any possible theory, i.e. scalar 
potential, that admits asymptotically AdS solutions. Instead we are interested in are 
theories that not only admit asymptotically AdS solutions, but also can be holograph- 
ically renormalized with local covariant boundary terms. The latter is a necessary 
condition for the locality of the holographically dual CFT and constrains the form of 




(4.45) 



e = e(A) 



(4.46) 
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the scalar potential. From the point of view of the algorithm presented in this paper, 
the requirement that the dual CFT can be holographically renormalized restricts the 
asymptotic form of the input function /i(p). 

The asymptotic expansions are qualitatively different for the cases e = 0, e = 1 
and < e < 1. Starting with the latter case one finds that the most general 
asymptotic form of the input function that appropriately generalizes the expansion 
(4.13) for generic < £ < 1 is 

n 

fi = 2p+A2n^-'p' + J2 A^n"^ p^""^ + Bn {log pY + Bn \ogp + cn + --- 

m=3 

(5.1) 

where the constants e, Ak, B, B, Q and C are the parameters specifying the input 
function. It should be stressed that even though these parameters are the input and 
can in principle be prescribed at will, this in fact does not always lead to a scalar 
potential of the desired form. In general there are some constraints that these input 
parameters must satisfy, which we will discuss momentarily. For < £ < 1 the 
parameter e defines a unique integer n > 2 via the inequalities 

2 

<!-£<-, (5.2) 



n + i n 

which tell us how many terms have to be included in the asymptotic expansion (5.1). 
The infinite set of special values of the parameter e, corresponding to the case when 
the first inequality in (5.2) is saturated, defines the series of so called resonant scalars 
(see e.g. [19] and references therein). 

resonant scalars: = n (5.3) 

For < e < I (or < e < ^, together with assuming vanishing A-s) the sum in the 
first line of (5.1) disappears and we are back to the case discussed in the previous 
section. 

Given the asymptotic expansion for the input function /i(p) we can obtain the 
asymptotic form of the function /2(p) by integrating (2.11) and that of the scalar 
0(p) by integrating (2.12). For /2 we get 

f2 = fi-n^ + --- (5.4) 

where j = QC, is again the angular momentum defined in (2.11). For the scalar field 
one finds 

0(p) = p^ + P30-P^ + ■ ■ ■ + Pn0"~V'" "2^' + /in+10"P"'logp 

+ <P+p-' + ■■■), (5.5) 
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where the coefficients (^_, /i^, /in+i and are determined in terms of the parameters 
of the input function as 



2B 



Z^ra+l 



= 



3(3e 



-03 



_ 2(l-2e) l-£ 2 2 



1 



aA_A I 



3(1 -e 

A I 



with 



a 



2A2e 



A, 



1 



ak 



Ak 



a" 



(5.6) 



(5.7) 



introduced for convenience. It is clear from these expressions that, as alluded to 
above, the parameters specifying the input function /i(p) cannot in fact be chosen 
completely at will. In particular, the inequality (2.13) requires A2 > 0.^ Yet another 
restriction concerns the parameter B. This parameter can be nonzero only in the case 
of resonant scalars since in all other cases it leads to a non-renormalizable theory. 

For e = the mass (4.9) saturates the BF bound and the asymptotic expansion 
of the input function /i(p) takes the form 



fi = 2p + An([og pf + sn(iog pY + Bn\ogp + cn + 

Integrating (2.12) in this case leads to 

0(P) =P"'/' (logp</)_ + </)+ + ■•■), 

where 



(5.^ 



(5.9) 




Finally, for e = 1 (vanishing mass) the input function /i (p) must be of the form 



/i = 2p + cn 



p 



(5.11) 



leading to 



= 0_+0+p-^ + --- (5.12) 
where 0+ = yJ—DVt and 0_ is the arbitrary integration constant of (2.12). 

^Thc case = is special since the asymptotic expansion for the scalar is apparently ill defined 
in this case. However, one should keep in mind that the above expressions for the coefficients of the 
asymptotic expansions are in general expressions for A^, B and B as functions of the parameter (/)_ 
and they remain valid even for = 0. For nonzero or A2, we can invert these relations and 
treat Ak, B and B as input, but this is not possible when A2 = 0. We conclude that if A2 = 0, then 
all Ak>2, B and B must be zero in order to obtain a scalar potential that admits an AdS vacuum, 
while </)+ is arbitrary. This case implies Brown-Henneaux boundary conditions = 0. 
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5.2 Scalar potential 

The asymptotic expansions of the previous subsection were constructed by demand- 
ing that the corresponding scalar potential V takes a certain asymptotic form so that 
the theory can be holographically renormalized. In this subsection we determine the 
potential V that follows from the asymptotic expansions above. 

Again, one has to distinguish the cases e = 0, e = 1 and < e < 1. In the latter 
case, the most general scalar potential that follows from (5.1) takes the form 

where the mass is given by (4.9) and the next couple of terms are 

I/O .^2 2,^ 3(3e-l)2(17£-7) ^ {le-2){l-ef 

8^ ; ^ 4 gv ; 4 128(1-6) ^ 245 

(5.14) 

Note that the cubic coefficient vs is non- vanishing only if the expansion coefficient A3 
is non-zero. By contrast, the quartic coefficient is in general non-zero even when 
both and A4 vanish. Moreover, at order n + 1 there is a difference depending 
on whether the inequality (5.2) is saturated or not. Namely, in the case of resonant 
scalars the coefficient Vn+i contains an additional contribution relative to the non- 
resonant case, given by 

-2e{l - e)fin+i. (5.15) 

By contrast, in the case e = the scalar potential is much less constrained and 
takes the generic form 

n0) = -^0' + o(0') (5.16) 
while for 6=1 the potential vanishes, leading to the case discussed in section 2.4. 

5.3 Holographic renormalization 

Now that we specified the scalar potentials we are interested in, we determine sys- 
tematically the counterterm function U{(j)) in (4.22) in this subsection. We can do 
this directly at the level of the minisuperspace model (2.4) using a Hamiltonian lan- 
guage where the radial coordinate plays the role of Hamiltonian 'time'. In the full 
Einstein-scalar theory this amounts to an ADM [20] formulation of the dynamics in 
the radial coordinate. 

From (2.4) (with S = j dp L) we obtain the canonical momenta 

f)r fjj 

n, = = Vol e-\jP = 771 = -Vol . (5.17) 

aA* ocp 

At the same time, these momenta are related to the on-shell action, or Hamilton's 
principal function S[X, 0], by 

n.^§ n^^f (5,18) 
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Combining these two expressions for the canonical momenta allows us to write the 
constraints (2.5c) and (2.6) as PDEs for the on-shell action S, namely 




2 + y(0) = 



(5.19) 



(5.20) 



Hamilton- Jacobi theory is based on the fact that finding a complete integral (i.e. a 
solution that contains as many integration constants as there are fields) of these 
PDEs is equivalent to completely solving the EOM. Mathematically, this follows 
from the theory of first order PDEs and the fact that Hamilton's equations are the 
characteristic equations of the Hamilton- Jacobi PDEs. 

The boundary term that is required to render the on-shell action finite agrees 
with a suitable solution of the Hamilton-Jacobi equations up to terms that vanish 
at infinity [21]. In fact, this same boundary term is required to make the variational 
problem well posed [22,23]. As we shall see, not any solution of the Hamilton-Jacobi 
equations can be used in the boundary counterterms, but the suitable solution is not 
unique either in general. The non-uniqueness originates in the freedom of choosing 
the integration constants in the solution of the Hamilton-Jacobi equations. These 
integration constants only affect the finite part of Hamilton's principal function. 
The holographic interpretation of this ambiguity is the usual renormalization scheme 
dependence of the renormalized generating functional in the dual field theory. This 
observation is particularly important for the problem at hand, because the angular 
momentum Jj corresponds to one of these integration constants and so it only affects 
the finite part of the on-shell action. For the purposes of determining the boundary 
counterterms, therefore, it suffices to consider only zero angular momentum solutions 
of the Hamilton-Jacobi equations. Any such solution takes the form 



The boundary term in (4.22) is exactly of the form (5.21), but with W{(j)) replaced 
with the function U{(j)). Note that W{(j)) and U{(j)) are not the same function, and 
this is why we use different symbols. U{(j)) is the function that defines the boundary 
counterterms, while W^cj)) denotes any exact solution of (5.22), i.e. the so called 'fake 
superpotential'. U{(f)) is in fact related to a particular solution W^cj)) in a way we 
will describe below. 

Before we analyze the various solutions of (5.22) for the scalar potentials we 
presented in the previous section, it is important to realize that the solutions W{(j)) 




(5.21) 



where W{(j)) is a function that is related to the potential by 



y(0) = 2 + 2W'^ - 2W\ 



(5.22) 
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of this equation determine the asymptotic form of the fields themselves via the first 
order equations 

X'^^Mx' ^=-Hipi (5.23) 

which are obtained by combining the expressions (5.17) and (5.18) for the canonical 
momenta. Since the point particle angular momentum Ji corresponds to a normal- 
izable mode of the metric, the asymptotic expansions obtained through (5.23) will 
be correct up to the order of this mode. These flow equations, therefore, apart from 
helping in determining which solution of (5.22) should be used in the counterterms, 
also relate the coefficients in the Taylor expansion of the potential to the coefficients 
of the asymptotic expansion of the input function /i(p). 

For < e < 1, the solutions of (5.22) with the potential (5.13) take one of the 
two possible asymptotic forms 

r 1 + %L02 ^ . . . 

[ 1 + ^ ^30^ + W40^ + ■ ■ ■ + + w)„,+i0"+Mog + A0^~ + ■ ■ ■ 

where the coefficients Wk are related to the coefficients in the potential, except 
for the constant A which is an integration constant of the first order equation (5.22). 
The dots in (5.24) stand for subleading terms that do not contribute to the on-shell 
action once the radial regulator is removed. Moreover, the logarithmic term can be 
present only in the case of resonant scalars (5.3). The function f/(0) that defines the 
counterterms must, therefore, take one of these two forms, with the terms denoted 
by dots dropped since they do not contribute to the on-shell action. To determine 
which of the two asymptotic forms f/(0) must take, notice that the coefficient of 
0^ in W{(f)) determines, via the first order equations (5.23), the leading asymptotic 
behavior of the corresponding solutions for 0. Namely, the two possible signs lead 
respectively to 

^ p-^±/2. (5.25) 

These are precisely the asymptotic forms of the two possible independent solutions 
for the scalar field in (4.16) and (5.5). The plus sign, therefore, corresponds to 
solutions of the EOMs with 0_ = 0, i.e. Brown-Henneaux boundary conditions. 
Hence, in order to ensure that the boundary term makes the on-shell action finite 
on all possible solutions the function ?7(0) must be of the second form in (5.24) [24]. 
Moreover, locality of the boundary term means that the term proportional to A can 
be included only when the exponent 2/A_ is integer, i.e. only for resonant scalars 
(5.3). In that case, A can be non-zero and its value can be chosen at will as this 
term only changes the finite value of the renormalized on-shell action and of the 
other renormalized quantities. It corresponds to part of the renormalization scheme 
dependence of the dual theory. There is a unique value. A*, of A in f/ (0) for which the 
renormalized stress tensor in the case of resonant scalars remains of the form (4.28), 
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and we assume that this choice of scheme is made. Finally, the logarithmic term 
in the case of resonant scalars would seem to make the boundary term non-local in 
0. This is a typical signature of a conformal anomaly. Locality can be preserved at 
the expense of introducing cut-off dependence in the boundary counterterms. Since 
~ p"^"/^, putting everything together, we conclude that for < £ < 1 the function 
U{(f)) in (4.22) takes the from 

f/(0) = l + ^02^W30'+u;40' + ---+u;„0"-^^w)n+i0"+'logpoo + A*0^ (5.26) 

where poo is the radial cut-off and the last two terms are present only in the case 
of resonant scalars. The coefficients Wk and Wn+i can be related via the first order 
equations (5.23) to the coefficients of the asymptotic expansion for the input function 
/i(p). For the first couple of terms we find 

1 , , 1, , 15(3e-l)2 ^ (1-ef 

while 

Wn+i = -ejin+i- (5.28) 

Note that the counterterms (5.26) are local in the induced fields, which ensures that 
the dual CFT is a local renormalizable field theory. It was precisely this requirement 
that led us to the form (5.13) of the scalar potential. 

For e = 0, the potential (5.16) leads to the two asymptotic solutions of (5.22) [18] 



a2 



4- y^ + J^*)+'' 

Again, the first solution corresponds to the mode (/>_ in (5.9) being switched off and, 
hence, the function U{(j)) must correspond to a solution of the second type in (5.29). 
As in the case of the resonant scalars, in order to preserve locality we must introduce 
explicit cut-off dependence in the boundary counterterms. Since log cj) ~ — ^logp, 
we deduce that the counterterm function for e = takes the form [24, 25] 

U{cP) = l + ]4>^(l--^), (5.30) 



4 V log po 

where poo is the location of the cut-off. 

Finally, for e = 1 {V{(j)) = 0) the scalar does not contribute to the divergences of 
the on-shell action and one only needs to remove the volume divergence by choosing 
in (5.22) the solution 

U{(j)) = l. (5.31) 



-34- 



6. Black hole thermodynamics 



Black hole thermodynamics provides some key insights into semi-classical and quan- 
tum gravity. Of particular importance is the black hole entropy and its microscopic 
description, for instance from a CFT perspective by virtue of the Cardy formula. In 
this section we address these issues for EDG. 

In section 6.1 we discuss basic thermodynamical quantities, like temperature, 
entropy, etc. We derive a formula for entropy in terms of curvature invariants. In 
section 6.2 we focus on black hole families and their associated solitons. We prove 
that these solitons can never have a (free) energy lower than global AdS for marginal 
boundary conditions. In section 6.3 we show the validity of the Cardy formula in 
EDG, again for marginal boundary conditions. 

6.1 Basic thermodynamical quantities 

Black hole thermodynamics can be studied efficiently on the gravity side in the Eu- 
clidean path integral approach, by exploiting results from the previous sections. We 
refer to [22] for a general analysis of black hole thermodynamics in the context of 
holography in asymptotically locally AdS backgrounds. This section is a particular 
application of the results there to the 3-dimensional case we are discussing here, tak- 
ing into account the modifications required by the possibility of generalized boundary 
conditions. See also [26] for a similar analysis of black hole thermodynamics in 2- 
dimensional EDG. 

Temperature Demanding the absence of a conical defect at the horizon p = ph 
leads to a periodicity in Euclidean time t ^ r + (3 that is identified with the inverse 
temperature. From the line-element (2.28) we obtain 



/3 = T- 



RdpN^ 



2nRh , , 

(6.1) 



Ph 3 

Notice that for black hole solutions the particle angular momentum j is positive. 

Entropy To derive entropy from scratch we could evaluate first the on-shell action, 
extract from there the free energy and then obtain entropy by taking the appropriate 
partial derivative of the free energy with respect to temperature. We shall provide 
these results below. Alternatively, for a two-derivative gravity theory it is known 
that the black hole entropy simply is given by the Bekenstein-Hawking formula 

S = AirAh (6.2) 

where Ah = 2TxRh is the horizon area. Using (6.1), an interesting reformulation of the 
result for entropy is 

S = 2Yo\Ei (6.3) 
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where Vol£; = J drdip. The entropy formula (6.3) highhghts the physical importance 
of the constant of motion j from (2.11). 

Our result (6.3) allows a relation between entropy and the ratio of curvature 
invariants, the square of the Cotton tensor and the square of the tracefree Ricci 
tensor, that resembles Penrose's Weyl curvature conjecture in four dimensions [27]: 



The prefactor is given by ^ = / d^XA/— 7, where the integral goes over some 
p = const, hypersurface with induced metric 7. The relation (6.4) is independent 
from the particular choice of hypersurface and can be checked easily using (2.19), 
(2.20) and (6.3). 

Angular velocity Since a black hole is a localized object, its angular velocity is 
measured with respect to that at infinity, which is a property of the asymptotic 
region 

uj = N^h-N^^- (6.5) 

In the AdS case we are currently considering and without loss of generality, we then 
fix the asymptotic form of the spacetime by requiring the boundary metric to be 
Minkowski. Namely, N^oo = and R^ = 2. This reduces the map between the local 
and physical coordinates (2.15) to (2.47), as was already used for the derivation of 
BTZ solutions. This restriction yields 

uj = tanhr; (6.6) 

with rj as defined in (2.47). 

Mass and angular momentum The holographic conserved charges associated 
with a boundary conformal Killing vector /c* are given by 

Q[k] = j dip {Tu) . (6.7) 


Since we know the renormalized holographic stress tensor, it is straightforward to 
evaluate these conserved charges. The black hole mass is obtained for k = dt 

2tt 



M = j d(f ((Too) cosh^ r] + (Tn) sinh^ r]) {6.1 




where the stress tensor should be evaluated with the appropriate boundary condi- 
tions, and the subindices here refer to the local coordinates (x°,x^). Substituting 
(4.28) for Dirichlet, or (4.37) for Neumann or Mixed boundary conditions, we obtain 

M = 271 j cosh2r/ - 2n(ri) . (6.9) 
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The black hole angular momentum is obtained for fc = 9^. 

27r 

J = j (Tt^) = 71 ((Too) + (Til)) sinh277 = 27ij smh2rj (6.10) 



Gibbs free energy and the first law of black hole mechanics The renormal- 
ized Euclidean on-shell action gives the Gibbs free energy [28] 

I = l3G{T,uj) (6.11) 

where / = —T is the renormalized Euclidean on-shell action and 

G{T,uj) = M -TS -ioJ (6.12) 

is the Gibbs free energy. This coincides with the Helmholtz free energy for vanishing 
angular velocity and angular momentum. Evaluating the Euclidean on-shell action 
we get 

G{T, to) = -271 j - 271 {Ti) . (6.13) 

Indeed, using the relations derived above, this expression reproduces the right hand 
side of (6.12), and —dG/dTl^^ reproduces the Bekenstein-Hawking entropy (6.2). 

Finally, as is shown in general in [22] these quantities satisfy the first law of black 
hole mechanics 

dM = T6S + uj6J (6.14) 

where the variations are taken at fixed source and fixed boundary condition for 
the scalar field. 

6.2 Black hole families and their solitons 

We have seen in previous sections how to systematically construct black hole families 
by varying, besides the rotation parameter t], the integration constant Q. We analyze 
here their thermodynamic properties. 
The black hole temperature, (6.1), is 

(6.15) 




27T cosh ?7 V 

with the value of Q for the unique static soliton solution associated to each 
black hole family, given in (2.46). As we have discussed, only marginal boundary 
conditions with vanishing source allow to freely vary the parameter Q. Namely, this 
is the boundary condition for which the {Q, ?7}-black hole family can be interpreted 
as thermal states in the same boundary CFT. Since for marginal boundary conditions 
the trace of the stress tensor is proportional to the sources, (4.40), the field theory 
energy and angular momentum are given by 

M = 87rn cosh 2r] J = Stt^ sinh 2r] (6.16) 
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We have chosen fl satisfying (2.44), such that j = 4f2 for simphcity. The entropy is 

5 = 327r2coshr/ v^^- (6.17) 

The charges of the sohtonic solution are 

Ms = -8nQs J = (6.18) 

Any of these 2-parameter famihes of stationary hairy black holes includes a 1- 
parameter family of extremal solutions. It is obtained by performing the double 
scaling limit ^2— )-0, |?7| — j-oo while keeping fixed ne^''''^0. We then have 

M=\J\=AttCi S = lQn'^^/nJl T = 0. (6.19) 

Stability of the hairy sector In addition to solutions with non-trivial scalar 
profile, the boundary condition (4.20) clearly allows solutions with vanishing scalar: 
the BTZ black hole family and its associated soliton, global AdS. The preferred 
solution will be that with smaller free energy for a given temperature. 

The free energy for hairy or BTZ black holes and their respective solitons is given 

by 

G = -87rQ Gs = Ms = -SnQs- (6.20) 

A first interesting implication is that inside each family there is a Hawking-Page 
phase transition at T=^. For T or equivalently fl<fls, the soliton geometry 
with a thermal circle is preferred, while for T > ^ the preferred one is the black 
hole. Comparing members of the hairy and BTZ families with the same temperature 
(6.15), we obtain 

G=^ Gbtz Ms = ^ M^,s (6.21) 

"AdS "AdS 

Which family dominates just depends on the quotient flg/flAds- 

We show now that Qg < fi^ds for any smooth solution with a center in EDG with 
non-trivial scalar profile. Relation (2.34) can be rewritten as 

fi, = i (6.22) 
2/ic 

where the subindex "c" indicates that the function should be evaluated at the center. 
We are assuming fi to be asymptotically positive, in particular /i = 2 at infinity. The 
inequality (2.13) requires then its second derivative to be negative for all values of 
the radial coordinate. Hence /i is a decreasing function, and we obtain 

fi, < ^ . (6.23) 

The inequality is saturated only when fi vanishes for all radial values, which implies 
that the solution is global AdS. We can state the following interesting result: any 
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smooth solution of EDG with vanishing trace of the boundary stress tensor has bigger 
mass than global AdS. 

Given that the free energy (6.20) is negative, the previous result ensures that 
a BTZ black hole has smaller free energy than the corresponding hairy one (for 
marginal boundary conditions). In spite of that, these hairy black holes are ther- 
modynamically perturbatively stable since their specific heat at constant angular 
potential is always positive: 



dT 



S > (6.24) 



This further implies that logarithmic corrections to the entropy from thermal fluc- 
tuations in the large VL limit take the form S'tot = 5 + 1 In S* + . . . , which is the same 
result as for the BTZ black hole, see e.g. [29] and references therein. 

We would like to stress, however, that the stability of a solution depends strongly 
on the boundary conditions. The same solitonic solutions may appear in other the- 
ories with different boundary conditions. The soliton mass Ms depends on them 
through a contribution from the trace of the boundary stress tensor, see (6.9). If 
that boundary conditions admit also global AdS as a solution, the soliton mass 
might be lower than that of AdS. Due to the inequality 

Ms = -8nQs - 271 (Ji) > M^ds - 27r (T^) (6.25) 

these solitonic solutions necessarily must have a positive trace of the boundary stress 
tensor, {Tl) > 0. 

6.3 Cardy formula 

The Cardy formula for specific hairy black holes in EDG was considered in [17]. 



5ca.dy = 47ry'-A+A+ + 47ry'-AoA- (6.26) 

where A^ = (M ± J)/2 are the Virasoro zero mode eigenvalues of the black hole 
states, and Aq are the lowest eigenvalues. The non-trivial question is what those 
lowest eigenvalues should be. Rewriting the entropy formula (6.17) in a suggestive 

way, 

S = 2ti^-Ms{M + J) + 2tt^/-Ms{M - J) , (6.27) 

we see that the Cardy formula is satisfied, S = Scardy, provided we identify Aq = ^ M^, 
as proposed in [17]. Since our conclusions hold for any (f2, ?7)-family of stationary 
black hole solutions we have generalized the results by Correa, Martinez and Troncoso 
from a set of specific examples to generic EDG. 
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7. Examples 



In this section we present some examples to illustrate our algorithm. In section 
7.1 we outline (and give explicit examples) how to obtain specific solutions that 
appeared in the literature on 3-dimensional gravity with a self-interacting scalar 
field [16,17,30-34]. In section 7.2 we present solitonic solutions that obey Brown- 
Henneaux boundary conditions. In section 7.3 we discuss a class of asymptotically 
AdS examples that violate Brown-Henneaux boundary conditions, for any scalar field 
mass in the window — 1 < < 0. In section 7.4 we consider solutions with expo- 
nential potential, which generally leads to solutions that are neither asymptotically 
fiat nor asymptotically (A)dS. 

7.1 Recovering known solutions 

By use of our algorithm we should be able to construct all stationary axi-symmetric 
solutions. In particular it should reproduce all solutions in the literature, most of 
which are static ones. To check this, given a static solution 

ds'^ = gtt dt"^ + g^^ d(f^ + Qrr dr"^ (7.1) 

there exists a coordinate p such that the metric looks like our (2.17). Taking /2 = 
—gtt, fi = Qipifi it follows that the required coordinate change is 

dr 1 

9ir) = ^ = , • (7.2) 

The EOM were written in the p coordinates. The explicit expression for r(p) might 
be hard to obtain, therefore we convert the relevant p-coordinate EOM into r coor- 
dinates (primes denotes derivatives with respect to r). 



(7.3) 



A V h 

Under the same change of coordinates, the integration constant j is 

2j=^?(r)(/^(r)/i(r)-/((r)/2(r)). (7.4) 

As examples we discuss in detail how to recover the solutions from [30] and [32] . 
The metric is given by 

_ {H + Bf _ 2 

" ~{H + By^^^^ ' ^'■^ " {H + 2ByF{r) ' " ^ ' ^^'^^ 

where 

H{r)=^-{r + ^MT4B^). ^ = ^ " + + ^) ' ^^'^^ 
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Using the coordinate transformation (7.2), 



AB + r 



9[r) = 

the scalar field is then determined as 

(b = A arctanhi 



B 



H(r) + B 



Taking into account the different normalizations, 4$ 
field used in [30], we arrive at the potential 

1 



(7.7) 



(7i 



!), where $ denotes the scalar 



( cosh*" $ + 1/ sinh^ $) 



(7.9) 



Furthermore, we find the relation between constants of motion and parameters 



J 



3B'(1 



(7.10) 



Since u is part of the potential and B^ is not, the quantity B^ plays the role of our 
constant of motion j, whereas is a free parameter of the model. 

Now we check that our general results for physical observables like mass and 
entropy reproduce the results obtained in [30] through the formalism of asymptotic 
symmetries. To this end we write the action of [30] in a way in which the different 
normalizations are more apparent: 



S 



1 

1 



d X^ 



-9 



R 
16 



where we defined 



16nG 



V{<P) 



2 

i? + 2 



(V$)2 
1 



47 + 15z/ 



(7.11) 
(7.12) 

(7.13) 



8" 32" 30720 " 

which corresponds to e = ^ (resonant scalars) according to (4.8) and (4.9). Note 
that the two leading terms show a universal behavior expected from (5.14) for s = \ 
and 03 = 0. As explained in detail in section 4 observables of interest depend on 
the chosen boundary conditions (4.33). In order to calculate the mass we need to 
write 0(r) as in (4.16) and then identify a function >V'((/)_) such that the regularity 
conditions (4.44) holds. In the case at hand e = 1/2 and therefore A_ = 1/2 and 
A+ = 3/2. 

Since we are only interested in the asymptotics it is sufficient to take the leading 
terms of p = p(r). Expanding and integrating (7.7) leads to 



- 2Br + 



P- 



(7.14) 
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Inverting the above expression and substituting into (4.16) gives 



0(p) = 2^/^v^p-V4 _ 1^21/^53/2; 



-3/4 



(7.15) 



so that we can identify 
identify 



1^1^^ and 



k _ 10 

48^" 



2i/4^3/2_ Using (4.44) we further 



(7.16) 



10^2 _|_ rpj^g integration 



or W(0-) = - Ja*^- + ^- Equivalently yV(0_) = 3 
constant w corresponds to an unphysical additive constant of the effective poten- 
tial (4.45). Assuming Dirichlet boundary conditions the renormahzed stress energy 
tensor (4.28) is 

(Too)+ = ^(9(l + z/) + 20). (7.17) 



For mixed boundary conditions it follows from (4.37) that 



(Too)- 



3B'(l + iy) + 2w. 



(7.18) 



To obtain the mass we integrate in ip (multiply by 27r) and use the reference's different 
normalization (divide by 167r), and obtain 



M =Ib\1 + u) + ]w 
8 4 



(7.19) 



which agrees with the result of [30] for w = 0. 

There exits a second class of solutions obtained from the potential (7.9) for z/ < 
presented in [32]. The metric is given by 



9tt = -r , 

Qrr 

fir) 



p(P-B{B-2f)^-B^u) 



(7.20) 
(7.21) 

(7.22) 



r-B( 



r - 5(v/^- 1))2 + 45v/^r) . (7.23) 



We again transform the coordinates via (7.2). The scalar field asymptotically is given 
by 



'R /R\3/2 



) 



and agrees with 



$fr) = arctanhi 



(7.24) 



(7.25) 
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by taking into account the different conventions, 4$ = (p. Note that in this case 
j = and we have = —3/4, = 3/2 and A_ = 1/2 as before. The change of 
variables from r to p is (perturbatively) 

= ^ + + - ' (7-26) 



r = 25 + V2V252 + P, (7.27) 

(f){p) = 2^/^V^p-V4 _ 1(21/453/2(7 _^ 3^))p-3/4 ^ (7 28) 

O 

leading to 

0+ = -1(2^/^5=^/2(7 + 3v/^)) , (7.29) 

0_ = 2^/^v/5 . (7.30) 

Through the regularity condition (4.44) we obtain 

W(0_) = --^(7 + 3v/^)0i^ (7.31) 



384 

Inserting everything into (4.37) gives 

2 



96 

W(0_) = --l-(7 + 3v/^)0i +w. (7.32) 



(Too)+ = ^B\7 + 3V^) (7.33) 
(Too)_ = (Too)+ + 2W(0_) = 2w (7.34) 

in agreement with [32]. Because the stress tensor does not depend on j (or equiva- 
lently B), this is called a degenerate state. 

7.2 Brown— Henneaux solitons 

Here we focus on asymptotically AdS solutions with specifically designed properties. 
We are particularly interested in regular solutions without black hole horizons. There 
are two possibilities to achieve this: either we make sure there is no horizon but rather 
a center, or we ensure that the horizon corresponds to a Poincare patch horizon and 
not a black hole horizon. We start with the latter option. 

Poincare patch soliton As discussed in section 2.2 static solutions exhibit a 
Poincare patch horizon if the functions fi and /2 simultaneously have a zero. This 
can be achieved most easily if both functions coincide with each other. A simple 
representative choice is 

^' = = 2" (l + ^) • 
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Figure 2: Potential 1/(0) for (7.37) 



with some positive fi. This choice ensures asymptotic AdS behavior with Brown- 
Henneaux boundary conditions for large p and a Poincare patch horizon at p = 0. 
The scalar field (2.12) is given by 



= TT -4arctan^— (7.36) 

where we chose the integration constant such that limp_!.oo = 0. At the Poincare 
patch horizon the scalar field approaches the value limp^o0 = — vr/-\/2. Evaluating 
the scalar potential (2.14) yields a simple result. 

^(^20) = -6sinV (7.37) 

Figure 2 depicts the potential (7.37) in the whole range of 0. The asymptotic 
AdS region corresponds to the maximum of the potential in the center, while the 
minimum corresponds to the Poincare patch horizon. The solutions following from 
the Ansatz (7.35) are a special case (j = 0) of the next family of solutions. 

Soliton with center As discussed in section 2.2 static solutions exhibit a regular 
center if the function /i has a zero, but /2 has no zero (or vice versa if the meaning 
of time and angular coordinate are exchanged). Let us start by choosing again 



P 

Then we can integrate the ODE (2.11) and obtain 



For ^ = the Poincare patch soliton of the previous paragraph is obtained. If ^ is 
positive the function /2 has a zero at some positive value of p. If C, is negative the 
function /2 has no zero in the whole domain p G [0, oo) and regular solitons with 
center are obtained. 
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The scalar field is given by (7.36). Evaluating the scalar potential (2.14) general- 
izes the result (7.37) to something more complicated. Consistently with our general 
discussion, the potential (7.37) is independent from the constant of motion Q, but 
like in the asymptotically flat example in section 3 it does depend on ^. 

It is straightforward to generalize the discussion above to functions with similar 
behavior, for instance fi = 2p (l + r2/(p+f2a^)), with some parameter a. The general 
features remain the same, but the potential is more complicated, which is why we 
refrain from presenting further explicit examples of this type. 

7.3 Non-Brown— Henneaux solutions 

An interesting class of input functions /i that asymptotically behave as (5.1) for 
^ e < 1 is given by 

A(p)=iViv^Y^(^)'^ (7.40) 

with A'^i = — 2=^7r/r [jzr^) ■ The advantage of the choice above is that we can exploit 
properties of Bessel functions Y„ and J„. Integrating the ODE (2.11) leads to 

e-l 

/2(p) = /i(r)+eiV2V^J^(^) " (7.41) 

with N2 = -2T^r(^)/(l -e). 

This spacetime has infinitely many Killing horizons and a curvature singularity 
at p = 0. The scalar field is given by a remarkably simple expression, 

^(P)=(^)*, (7.42) 

which allows a straightforward inversion to extract p{(f)). The expression for the 
potential for different values of e is lengthy and not too illuminating, but can be 
obtained analytically. For example, if e = | and = we obtain 

V{<P) = 2 + ^ [2cj>\<f^' - 54) Yl{<P) - - 216) YMM^I^) 

- (30'^ - 17602 + 1728) Y2(0) + 2</.^Y^(0) - (p^Yl{(f))] (7.43) 

As expected, the potential does not depend on Q. Note that there is a subtlety 
in obtaining the full potential as a function of the scalar field for non-integer 
because the potential becomes complex for < 0. Remember that we fixed a sign 
ambiguity in (2.12) by assuming that is non-positive. By allowing also to be 
positive we can glue together the two branches and obtain the full (real) potential. 

In figure 3 we plot the potential for three different values of e and also show the 
oscillatory behavior in the deep infrared (large values of 0) which originates from the 
behavior of the Bessel functions. 
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7.4 Monomial Ansatz — solutions with exponential potential 

Here we study a simple case with V ^ that is neither asymptotically flat nor asymp- 
totically (A)dS, just to demonstrate that such solutions exist and can be constructed 
as well following our algorithm. They are of little physical interest, though. 
Consider the power-law Ansatz 

fi = 2p{p/nr-' 0<a<l (7.44) 

with some Q 0. Then /2 (as long as a 7^ |) is determined as 

f2 = xfi + Y^(p/^y"' (7-45) 

where x is some constant. The scalar fleld is given by 

(j) = (j)o± ^/a{l-a) In ^ . (7.46) 

The potential is then given by 

V{(j); a, vo) = 2- vq e^'* . (7.47) 

Note that V is independent from the constants of motion j = and Q. Upon 
redeflning 0o also the constant x can be absorbed. The constants a and ^0(0;, (po) 
are then parameters of the model. No further constants of motion appear in this 
construction, which means that further solutions to a model with the potential (7.47) 
must have non- monomial solutions for fi. The case Oi = \ can be treated similarly, 
leading to the potential (7.47) with a = \ and to /2 = /i (x + | lii(p/^))- 

Non-trivial solutions with exponential potential always have a scalar fleld that 
diverges logarithmically at p — )■ 00 and p — t- 0. 
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8. Generalizations 



Our algorithm to bootstrap gravity solutions by providing geometric input, in the 
form of one free function /i(p) subject to the inequality (2.13), in principle allows to 
construct all stationary axi-symmetric solutions of all EDG models in three dimen- 
sions (2.1). Finding all solutions for fixed scalar potential V turned out to be very 
hard, so we restricted ourselves to a co-dimension 2 family of solutions. 

Here we discuss some possible generalizations of our methods. We start with gen- 
eralizations within 3-dimensional EDG. It is straightforward to drop axi-symmetry 
and replace it by translation invariance, which leaves the local analysis unchanged 
and just affects some of the global properties. For instance, the hairy black hole 
solutions in section 3 then lead to hairy black brane solutions. In section 2 we have 
repeatedly assumed certain smoothness conditions, in particular analyticity. It could 
be interesting to relax these assumptions and check if it leads to qualitatively new 
solutions. A precedent for this is the numerical evidence for stationary axi-symmetric 
solitonic solutions of topologically massive gravity that are not analytic [35]. The 
solitonic solutions discussed in sections 2.3 and 6.2 are special and can never have 
a mass below that of AdS for marginal boundary conditions. However, it is known 
from "designer gravity" [2] that solitons can exist whose mass is below that of AdS. 
These solitons are in fact equivalent to our solitons as solutions to the EOM, but 
they require non-marginal boundary conditions and must have a positive trace of the 
associated boundary stress tensor. 

Within EDG there are no locally Schrodinger or Lifshitz solutions, because in 

•"2 

terms of Clement's Ansatz (2.3) these solutions require 7^ X and X = 

[4,35]. There is no solution to the Einstein equations with these properties, since 
(2.5a) implies proportionality of X and X. This provides a motivation to consider 
theories that are more general than EDG. 

It would be of interest to consider Einstein-Maxwell theory along the lines of 
EDG, since this theory has the same number of physical degrees of freedom as EDG, 
so one might expect a comparable amount of complexity. The crucial differences to 
EDG are the absence of some free function (no scalar potential) and the presence 
of a gauge symmetry. For applications to holographic superconductors [36, 37] it is 
of interest to combine these two theories and study Einstein-Maxwell theory with 
a charged scalar field. This may allow to find analytic toy models for holographic 
superconductors. Adding Chern-Simons terms — either gravitational ones or for a 
Yang-Mills field [38] — provides another interesting generalization that is capable of 
producing locally Schrodinger, Lifshitz and warped solutions. 

In higher dimensions, D > 3, similar methods can be used, provided one has 
sufficiently many, D — 1, mutually commuting Killing vectors available. Then one 
can torically reduce the action, derive an analog of the Clement action (2.4) and 
build up an algorithm analog to the one described in section 2. 
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